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A general multi-orbital Hubbard model, which includes on-site inter-orbital electron hoppings, is 
introduced and studied. It is shown that the on-site inter-orbital single electron hopping is one of the 
most basic interactions. Two electron spin-flip and pair-hoppings are shown to be correlation effects 
of higher order than the on-site inter-orbital single hopping. It is shown how the double and higher 
hopping interactions can be well-defined for arbitrary systems. The two-orbital Hubbard model is 
studied numerically to demonstrate the influence of the single electron hopping effect, leading to a 
change of the shape of the bands and a shrinking of the difference between the two bands. Inclusion 
of the on-site inter-orbital hopping suppresses the so-called orbital-selective Mott transition. 

PACS numbers: 71.10.-w, 71.10.Fd, 71.27.-|-a, 71.30.-|-h 



Strongly correlated electron systems have drawn great 
interest, as they exhibit many exotic phenomena, such as 
metal-insulator transitions, unusual forms of m agn etism, 
superconductivity, and heavy-fermion behavior.lJJ^ These 
systems are therefore investigated theoretically in both 
model and ab initio calculations. One of the most es- 
sential tasks is to model the system appropriately, i.e., 
such that a realistic physical picture can be attained. 
In condensed matter theory the Hubbard model, which 
was originally proposed in the early sixties,'2HS! is one of 
the simplest, yet also the most important and most fre- 
quently studied lattice model to investigate strongly cor- 
related electron systems. It sets up a competition be- 
tween an inter-site quantum mechanical hopping term 
and an on-site Coulomb interaction term. As a conse- 
quence the model can describe various non-trivial phe- 
nomena. Due to its simplicity and because the model 
has captured the essence of strongly correlated electron 
systems, the Hubbard model has been widely used:- '^ 

In a realistic situation an atom in a correlated sys- 
tem will usually have several partial filled orbitals and 
should therefore be described with a multi-orbital (MO) 
Hubbard model. With the inclusion of orbital degrees 
of freedom, inter-orbital interactions have to be included 
for such a system. The inter-site hoppings should now 
be a sum of the inter-site intra-orbital hoppings and the 
inter-site inter-orbital hoppings, where the two kinds of 
hopping are defined according to the orbital indices of 
the start and the end orbitals for hopping electron are 
identical or not, respectively. Moreover, besides the re- 
main competition of the intcr-site hopping and on-site 
Coulomb interactions, one new kind of on-site interac- 
tions will join into this competition, i.e., the on-site 
inter-orbital hoppings. Acting as the inter-site hoppings, 
the on-site inter-orbital hoppings are also associated to 
the kinetic energy of electrons. All the inter-site hop- 
pings, on-site inter-orbital hoppings and on-site Coulomb 
interactions (which also divides as on-site intra-orbital 
Coulomb interactions and on-site inter-orbital Coulomb 
interactions) constitute a competition for various interac- 
tions existing in a realistic multi-orbital system. There- 



fore, for an arbitrary multi-orbital system, we propose a 
generalized MO-Hubbard model as follows, 

H = ~ ^ tjjhnfiicrfjma + Uunj^flu^ 
ijlma,i^j il 

ilmaa' ,l<rn 

+ ^ ] {-^ibnafima filer ^ -^ilmcr fii^fima) T i^) 
ilm(7,l<m 

where i,j are site indices, /,m are orbital indices, and 
a, a' are spin indices. The , /i/o- are the creation and 
annihilation operator, respectively, for spin <j in Z-th or- 
bital on site i. The first term is the inter-site hopping 
term, where tijima is the hopping amplitude for spin a 
hopping from m-th orbital on site j to Z-th orbital on site 
i. The second term is the intra-orbital Coulomb inter- 
action and the third term is the inter-orbital Coulomb 
interaction, where Virn-\^—Vim^-\ + J and J is the Hund's 
coupling constant. The last two terms in the Hamilto- 
nian are the on-site inter-orbital single electron direct 
hopping terms, where the and luma are the hop- 

ping amplitudes, as tijim in the first term. Note that 
this part is written in a form of one term and its conju- 
gate, so that one can easily observe that the system is 
in a statistical equilibrium state. For the half filled case, 
/limt = hi-m\, = l*ira^ = ^ira^^ the model automatically 
obeys spin rotational invariance. Here the on-site inter- 
orbital single electron hopping is to be considered as one 
basic interaction along with the intcr-site hopping and 
the on-site Coulomb interaction. One may argue that 
with specially chosen basis the on-site inter-orbital hop- 
ping will disappear. However, we are here discussing a 
general multi-orbital system with arbitrary basis. It is 
convenient to investigate such a system because in some 
experiments the natural chosen basis in a material may 
not be orthogonal. Thus studying a system theoretically 
with a same basis set as in the experimental material 
can easily gives the results comparative to the physical 
observables obtained in experiments. Moreover, we find 
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that the inclusion of on-site inter-orbital single hoppings 
has introduced some new and interesting theor y. 

Recently we have numerically studiecP^EI] the first 
three terms in Eq. (flj within the dynamical mean field 
theory (DMFT).!^ i^' in this work we formulate theoreti- 
cally and evaluate numerically the influence of the on-site 
inter-orbital single electron hopping. Using the DMFT 
the proposed generalized Hubbard model can be mapped 
to a generalized single impurity Anderson model (SIAM) 
along with a self-consistency condition. The Hamiltonian 
of this generalized SIAM is 

kla la I 

H~ ^ ^ Ulmaa' ^la^ma' 

Imaa' ,l<m 

+ X i^lkAkafl'y + Vlkaftcika) 
Ika 

+ ^ ] (^^Imafmafla^^lmiyfi^fma), (2) 
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where the first term is the energy of the conduction 
electrons (bath), c|j,^ and Qfeo- are correspondingly the 
creation and annihilation operators of conduction elec- 
trons. The second term is the energy of the localized 
electrons. The third (fourth) term represents the intra- 
orbital (inter-orbital) Coulomb interactions. The fifth 
summation is the hybridization term that gives the in- 
teraction between the bath and the localized electrons. 
The sixth summation represents just the on-site inter- 
orbital single electron direct hoppings on the impurity 
site, where for convenience we have dropped the site in- 
dices in the symbols. These inter-orbital single electron 
hoppings are on-site interactions which do not change in 
the mapping. We can hence study these terms in Eq. ([I]) 
by equivalently studying the mapped SIAM. 

We use the equation of motion (EOM) method to solve 
the DMFT impurity problem. In this method we com- 
pute the equations of motion according to the equation 

w<A;B>= ([A,B]+)+<[A,H.„p];B>, (3) 

where we have used the Fourier transform of the dou- 
ble time temperature-dependent retarded Green's func- 
tion (GF) < A{t');B{t) >, i.e., < > is defined 
in uj space.fi^ [, ]+ means the anti-commutator and [, ] 
the commutator. The first term on the right hand side 
(RHS) labels the interaction associated to <C A; i? 
the second term describes the involvement of the higher- 
order interactions, where the higher-order GFs will ap- 
pear. Calculating the EOM of these higher-order GFs, 
GFs of even higher-order will appear in the newly derived 
EDMs. Repeating this procedure, more higher and even 
higher order GFs will appear. Each GF is associated 
with an order and one physical interaction. The order of 
the GF approximately labels the weight of its associated 
interaction. Approximately, the GF of lower order asso- 



ciated with an interaction will give a larger contribution 
than higher-order GFs associated with the interaction. 

This statement can be demonstrated considering the 
following procedure. The 1*^* order EOM is 

(W + ^ - Efma) < fma, fLa > = 
l,l^m 

^ ^la' Ima 5 I ma ^ ) 

^ ^ ^vaka ^^^vaka ^ frna^^ 
k 

- X </i-;/L». (4) 

where the last term is generated by the on-site inter- 
orbital single electron hopping terms, and ^ is the chem- 
ical potential. Next, we calculate the second-order EOMs 
of those newly appeared higher-order GFs on the RHS of 
Eq. 0, e.g., the EOM of the GF is 

(W /i - Sfla) < fla] fLa > = 

{[fia, /Lit) + Uii <^hi,>fi,; /L> 

+ X {Ul'laa -^ni'afla-jLa^ 

+ Ui'i„'„ <^ni,^,Jiy, fl^^ > ) 

+ X ^fe'^ -^Clka-Jla^ 
k 

- X «/;'<^;/L»: (5) 

where we can observe from the derivative procedure that 
the last term actually reflects a physically consequent 

double hopping interaction, fl^{t2)fva{t2)f}na{h)fia{h)- 
Now calculating further the third-order EOMs of the 
GFs appearing on the RHS of Eq. ([5]), and taking <C 



fima' ficr fma ^ ^-^i iUustration we obtain 
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where V^{l,m) means Z't^Z and I'^m. The last six terms 
on the RHS are generated by the on-site intcr-orbital 
single electron hoppings. 
Next, noting that 

[fma, fma' fma flcr' flo-] ~ ~ flna' fl'^' fli^ ^ C^) 
[fma: f 1^1 fma fma' fla] f fma' fla ^ (8) 

we recognize that <- f^^, fia' fia] fLa'> in Eq. ^ is ac- 
tually associated with the pair-hopping term,— ""^ while 

fL' fma' fla', fL a ^ c orresponds to the so-called spin- 
flip exchange tcrmT^'^ One can note that both the pair- 
hopping and spin-flip exchange terms can be reproduced 
by on-site inter-orbital single hopping terms. The pair 
hopping is thus a special type of double hopping in which 
both the two electrons in the l-th orbital hop to the m-th 
orbital simultaneously. Similarly, the spin-flip exchange 
term is a certain kind of double hopping, in which one 
particle with spin a hops from the ^-th orbital to the 
m-th orbital as well as one particle with spin a' hops 
from m-th orbital to Z-th orbital at the same time, where 
(Ty^cr' . Thus we have given them a more general definition 
for arbitrary multi-orbital systems. Notably, in a multi- 
orbital system the possible double hoppings are not only 
the pair-hopping and the spin-flip exchange term. The 
other four terms in the last six terms on the RHS of 
Eq. ^ also correspond to certain forms of on-site inter- 
orbital double hoppings at an equal time limit. 

If we continue the procedure to calculate higher-order 
EOMs, more and higher-order GFs will appear which 
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FIG. 1: Illustration of possible simultaneous double hoppings 
in a multi-orbital system. The straight arrows label the spin 
up or down of an electron and the bend arrow gives the hop- 
ping direction, (a) the spin-flip exchange term; (b) the pair- 
hopping term; (c) two electrons in different orbitals with dif- 
ferent spins hop to the same orbital; (d) two electrons in one 
orbital with different spins hop to different orbitals; (e) two 
electrons in different orbitals hop to different orbitals. 



are associated with higher-order interactions. In view 
of the above, the on-site inter-orbital single hoppings 
fm'ai^'i)fma{ti) are the most basic physical interactions 
in Hubbard-like strongly correlated systems. These will 
furnish the higher-order double-time double-hopping in- 
teractions, fl„{t'^)fm:„,{t'i)fia{t2)fma'iti), that is. One 
single hopping occurs at time ti and another hopping 
occurs at time ^2- In the same manner higher-order hop- 
pings, as three-time three-hopping interactions and four- 
time four-hopping interactions etc. will appear. These 
multi-hopping interactions are embedded in our formula- 
tion and appear in higher-order EOMs. 

Next, we investigate in more detail the double-time 
double hopping at equal times, i.e., the two hoppings oc- 
cur simultaneously. The general form for the double hop- 
pings is fl^fm,„,fiafma', whcrc I' and m'^m. This 
equation gives all the double hoppings that possibly ex- 
ist in a multi-orbital system. When l' =m' l=m, cr'j^<J, 
it is the pair-hopping between orbitals.l- ^^ l ^^ l When l'=m, 
l=m' , a'j^cr^ it is the spin-flip exchange term. The whole 
set of on-site double hoppings include several kinds of 
interactions. If separated according to initial and final 
orbital of the hopping, we can schematically represent 
them as shown in Fig. [l] Writing 

fl'afm'a'fl'^fma' ~ ^ ^ fma fla' fl'^ fma' 

mlaa' .m^l,a^a' 

+ fla fla' fma fma' 

mlaa' .m^l,a^a' 

+ y ] {fmafma'fl'<^fl'^''^fl'afla'fmafma') 

mil' .1' ^l.a' ^a 

+ fl'afm'a'fl<^fma', (9) 

[ll'mm') all different 
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where the first term is just the spin-flip term and the 
second term is the pair-hopping term. Note that the 
fifth term can split into three sub-terms according to 
the two spins present in different orbitals. And, Fig. [T] 
only gives an illustration of the double hoppings, the 
terms should sum over all orbitals. If there are less or- 
bitals, the number of terms will be reduced accordingly. 
For three-hopping, four-hopping, and even higher-order 
multi-hopping interactions, a schematic representation 
similar to Fig. [T] can be made. 

From the above discussions, it can be recognized that 
the usually studied spin-flip and pair-hopping terms are 
higher-order consequences of the on-site inter-orbital sin- 
gle electron hoppings and Coulomb interactions. As the 
on-site inter-orbital single hopping is of lower order it 
will have a larger weight (or probability) than the spin- 
flip and pair-hopping terms. We note that, from the 
Eqs. (|4| and the off-diagonal GFs <-fia\fLa'> re- 
late to the on-site intcr-orbital single hopping. If one 
neglects the on-site inter-orbital single hopping but in- 
cludes the pair-hopping and spin-flip terms in a model 
Hamiltonian, the obtained off-diagonal GFs will be inex- 
act. Furthermore, in any self-consistent theory where the 
number of orbitals is larger than two, more inter-orbital 
double hopping terms have to be included, as shown by 
Eq. (|9]). Such terms can be well treated in our theory, 
because they will be automatically generated from the 
on-site inter-orbital single hoppings. 

The hopping amplitudes hma have the physical mean- 
ing of the hopping amplitudes. They are to be deter- 
mined for each studied system, but not fixed. For ex- 
ample, if in a certain system the a spin channel is fully 
occupied for one orbital all the incoming single hoppings 
of spin a into this orbital will be suppressed. In a param- 
agnetic system one can consider, for simplicity, that each 
orbital has the same incoming and out-going inter-orbital 
single hoppings. Importantly, all studies of higher-order 
hopping terms must be made on top of having the on-site 
inter-orbital single hoppings included. 

As an example we have numerically studied the two- 
orbital Hubbard model in paramagnetic case to show 
the influence of the on-site inter-orbital hopping. This 
interaction was explored in Ref. 1151 but not in a self- 
consistent manner. Theoretically, in the situation where 
the neighboring sites are identical to the impurity site, 
the influence of the on-site inter-orbital hopping is simi- 
lar to that of the inter-site inter-orbital hopping, as the 
orbitals on the neighboring sites are identical to the or- 
bitals on the impurity site. In order to cleanly distinguish 
the contribution of the on-site inter-orbital hopping from 
that of the inter-site inter-orbital hopping, we set the 
inter-site inter-orbital hopping and on-site inter-orbital 
fluctuations to zero. We used the MO-EOM impurity 
solver of Ref. 1101 in combination with genetic algorithm 
techniques .ESEH 

The computed quasiparticle densities of states (DOS) 
are shown in Figs. [2] and |3] for the cases that the two 
orbitals have identical or different band widths, respec- 
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FIG. 2: Quasiparticle DOS for the paramagnetic half-filled 
two-orbital Hubbard model on the Bethe lattice, computed 
for the case of identical band widths. The parameters used are 
the half bandwidth D2=D\ = 1, temperature r=0.01, and the 
Coulomb U as specified in the panels. Line 1) is obtained for 
the inter-site intra-orbital hopping parameters t2=t\=Di/2 
and no on-site inter-orbital hopping. Line 2) uses the same pa- 
rameters as 1) but the on-site inter-orbital hopping ti2=t\/2. 
Line 3) uses the parameter t^°^=t\ -1- ti2=Di/2 and ti=ti2. 
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FIG. 3: As Fig. [2) but for the case of different band widths. 
The left panels depict the DOS for the narrow orbital and the 
right panels that for the wide orbital. The parameters used 
are D2=2Di=2 and T^O.Ol. The hues 1) - 3) are obtained as 
explained in Fig. [2] only t2 is modified correspondingly with 
D2.— 'The narrow orbital's DOS for lines 2) and 3) with 17=3 
in (c) have been shifted vertically for visibility by 0.1 and 
0.05, respectively. 



lively. Line 1) is calculated without on-site inter-orbital 
hopping, line 2) with the same inter-site intra-orbital 
hopping and a nonzero on-site inter-orbital hopping, and 
line 3) with both inter-site intra-orbital and on-site inter- 
orbital hoppin g b ut the total hopping amplitude equals 
that for line l).'^Line 1) in Fig. ^ illustrates that there 
is a metal-insulator transition controlled by U. Line 2) 
shows that adding the on-site inter-orbital hopping in- 



5 



teraction the Hubbard bands are somewhat broadened. 
This happens because the on-site inter-orbital hopping 
effectively increases the total hopping amplitude so that 
the electrons gain some itineracy. In Fig. |3] line 1) shows 
that in the absence of on-site inter-orbital hoppings there 
IS an orbital selective Mott transition (OSMT) .EH How- 
ever, from both lines 2) and 3) one observes that the nar- 
row and wide orbitals simultaneously change from metal- 
lic states to insulating states along with the increase of U. 
Therefore, the OSMT shown with line 1) is suppressed 
with the inclusion of the on-site inter-orbital hopping. 

To summarize, we introduced a general MO-Hubbard 
model that can readily be employed to study a broad 
range of multi-orbital systems. We have shown analyt- 



ically and numerically the influence of introducing the 
on-site inter-orbital single hoppings. When these exist 
in a correlated electron system they will greatly change 
its properties. Higher-order effects, such as spin-flip ex- 
change and double hoppings, have to be studied on top 
of the on-site inter-orbital single hoppings, which, as out- 
lined, can be done in a well-defined way for arbitrary cor- 
related systems. The developed theory is expected to be 
beneficial for studies of unsolved correlation-related phe- 
nomena and to trigger inspiring theoretical studies and 
discoveries. 

We thank the Swedish Research Council (VR) and 
SKB for financial support. Computer time from the 
Swedish National Infrastructure for Computing (SNIC) 
is acknowledged. 



* Iqingguo .feng@physics.uu.se] 

^ N. Grewe and F. Steglich, Handbook on the Physics and 
Chemistry of Rare Earths, ed. K.A. Gschneidner, Jr. and 
L. Eyring (Elsevier, Amsterdam, 1991), Vol. 14, p. 343. 

^ M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

^ J. Hubbard, Proc. Roy. Soc. A 276, 238 (1963). 

" M. C. Gutzwiller, Phys. Rev. Lett. 10, 59 (1963). 

^ J. Kanamori, Prog. Theor. Phys. 30, 275 (1963). 

" J. Quintanilla and C. Hooley, Physics World 22, 32 (2009). 

^ H Tasaki, Prog. Theor. Phys. 99, 489 (1998). 

* T. Maier, M. Jarreh, Th. Pruschke, and J. Keller, Phys. 
Rev. Lett. 85, 1524 (2000). 

^ D. J. Scalapino, J. Supercond. Nov. Magn. 19, 195 (2006). 
^'^ Q. Feng and P. M. Oppeneer, J. Phys.: Condens. Matter 

23, 425601 (2011). 

Q. Feng and P. M. Oppeneer, J. Phys.: Condens. Matter 

24, 055603 (2012). 

W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 



(1989). 

A. Georges and G. Kothar, Phys. Rev. B 45, 6479 (1992). 
D. N. Zubarev, Sov. Phys. Usp. 3, 320 (1960). 
A. Koga, N. Kawakami, T. M. Rice, and M. Sigrist, Phys. 
Rev. B 72, 045128 (2005). 

P. Werner, E. GuU, and A. J. MiUis, Phys. Rev. B 79, 
115119 (2009). 

" X. Dai, G. Kothar, and Z. Fang, arXiv:cond-mat/0611075 
^® Q. Feng, Y.-Z. Zhang, and H. O. Jeschke, Phys. Rev. B 
79, 235112 (2009). 

Q. Feng, Ph.D. Thesis, University Frankfurt (2009). 
Specifically, the parameters used in Fig. [2] are: hue 1) 
ti=t2=Di/2, ti2=t2i=0; line 2) ti^t2^Di/2, fi2=t2i = 
_Di/4; line 3) ti=t2=-Di/4, ti2=t2i=-Di/4. The parame- 
ters used in Fig. [s] are: fine 1) ti=Di/2, t2=D2/2, ti2 = 
t2i=0; line 2) ti=Di/2, t2=I>2/2, ti2=t2i=Di/4; line 3) 
ti=Di/4, t2 = D2/2 - Di/4, ti2=f2i=£>i/4. 
21 A. Liebsch, Phys. Rev. Lett. 95, 116402 (2005). 



